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Abstract 

The Casimir surface force density Fona dielectric dilute spherical 
ball of radius a, surrounded by a vacuum, is calculated at zero tem- 
perature. We treat (n — 1) (n being the refractive index) as a small 
parameter. The dispersive properties of the material are taken into 
account by adopting a simple dispersion relation, involving a sharp 
high frequency cutoff at u = ujq. For a nondispersive medium there 
appears (after regularization) a finite, physical, force F nondlsp which is 
repulsive. By means of a uniform asymptotic expansion of the Riccati- 
Bessel functions we calculate F nondlsp U p to the fourth order in 1/u. 
For a dispersive medium the main part of the force F dlsp is also re- 
pulsive. The dominant term in F dlsp is proportional to (n — 1) 2 luq/ci , 
and will under usual physical conditions outweigh F nondlsp by several 
orders of magnitude. 
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1 Introduction 

Consider a compact dielectric nonmagnetic ball of radius a, surrounded by a 
vacuum. The permittivity of the material is e = n 2 , n being the refractive 
index. The purpose of the present work is to calculate the "ordinary" Casimir 
surface density F on the sphere (i.e. with the exception of the electrostrictive 
force), at zero temperature, making use of the same calculational technique 
as in an earlier work |IJ.. For mathematically simplifying reasons we will 
be concerned only with the case of dilute media, meaning that n - 1< 1. 
As in [|l| , we will make use of the uniform asymptotic expansion (also called 
the Debye expansion) for the modified Bessel functions [F2J. Here is 

the expansion parameter, with v = I + 1/2 , / = 1,2,3, In the main 

body of this paper we will expand the formalism to second order in {n — 1), 
and to fourth order in \jv . In Appendix A we carry out the expansion 
in [n — 1) three orders further, to the fifth order, at the price, however, 
of retaining the Debye expansion to second order in \jv . This treatment 
improves considerably on the theory given in [0]. We use the mathematical 
package Maple V Release 4 Power Edition. It turns out to be necessary 
to keep relatively high accuracy in this problem, since subsequent terms in 
the uniform asymptotic expansion alternate in sign. Thus, if one simply 
truncates the series after a few terms, without observing this alternating 
effect, one risks getting even the sign in the final force wrong. 

The ultimate goal in a calculation of this type is to get a Casimir force 
expression that is directly comparable to experiment; if not a real experiment, 
at least a Gedanken experiment. To achieve this goal, one must observe 
that there are several contributions to the surface force. Let us try to put 
the things into perspective, by writing down the general expression for the 
electromagnetic volume force density f in a charge - free nonmagnetic medium 
(cf., for instance, the book M or the review M) 



Here p is the mass density of the fluid. 

Consider the first term in (1). We may call it the Abraham - Minkowski 
force. It is equal to zero in the interior homogeneous region of the medium 
and acts in the boundary region around r = a only. ^From thermodynamic 
perturbation theory we know in general that e > 1 for a dielectric at zero 
frequency 0, Sect. 14, and we expect the same to be true for moderate and 
high frequencies also, except possibly from special frequency bands. That is, 
since the sphere is surrounded by a vacuum, we would expect the first term 
in (1) to give a repulsive surface force. This turns out to be true also, in 
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classical theory. In quantum theory, however, we cannot be sure in advance 
about the direction of the force; the result depends on the magnitude of the 
contact term being subtracted off in the regularization procedure. One has 
to carry out the calculation in detail in order to determine even the sign. 
The Abraham - Minkowski force can in turn be divided into two parts: 

(1) A nondispersive part. It corresponds to letting the ultraviolet (nondi- 
mensional) cutoff frequency xq go to infinity. The early works on the Casimir 
effect took into account this part of the force only. The nondispersive force 
is formally divergent, but can usually be easily regularized. The recipe is 
to (i) put Xq = oo ; (ii) evaluate the remaining divergent sum over I , from 
I = 1 to I = oo, by the Riemann zeta - function method. We stressed the 
efficiency of this recipe in ffl; cf. also some other related prior papers || - 
0. The method gives, in a simple way, results that are in agreement with 
what can be obtained in other ways. The quantum nondispersive force turns 
out in our calculation to be repulsive, thus qualitatively in agreement with 
the classical result noted above. 

(2) There is a dispersive part of the Abraham - Minkowski force, caused 
by the cutoff xq . In an earlier work on a spherical shell [[?[] we actually 
found the dispersive force to be attractive. The dispersive force is under 
usual physical conditions stronger than the nondispersive force. Candelas || 
seems to have been the first to emphasize the importance of the dispersive 
force. His general ideas were later on essentially supported by explicit model 
calculations ||, [TO ]. 

Consider next the second term in (1). This is the electro striction force. It 
generally acts inwards. Whether this force is detectable in a real experiment 
or not, depends on the detailed circumstances. Usually, there is established 
a compensating force from the elastic pressure in the interior of the fluid, so 
that the effect from electrostriction does not show up explicitly in the mea- 
surements. Only under special conditions, such as in the Goetz - Zahn non- 
equilibrium experiment testing the attractive force between two condenser 
plates immersed in a dielectric liquid when there is a high-frequency electric 
field between the plates [|TTJ, p!2| , does one have the possibility to measure 



the electrostriction force. The point is that the electromagnetic force will 
then vary so quickly that the elastic pressure does not have sufficient time to 
build itself up. What determines the time scale here, is obviously the transit 
time for sound waves in the body. (The Goetz - Zahn experiment, as well 
as some other related experiments, are discussed in detail in [|J.) In (1), we 
have deliberately written [dejdp)T as an isothermal partial derivative. This 
is often appropriate in practical cases, but is an assumption that has to be 
considered with some care since if the deformation of the body is occurring 
rapidly it will be more appropriate to take the adiabatic derivative {de/dp)s 
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instead. If the fluid is non- polar, as one can usually assume at high frequen- 
cies, then one need not distinguish between the isothermal and adiabatic 
cases since the permittivity depends on density only, though the Clausius - 
Mossotti relation. The electrostrictive contribution to the Casimir effect has 
not been paid much attention to in the past. There is a paper by one of us 
some years ago |13|], and recently Milton and Ng have returned to the topic 



14] . In the present case, where we are considering the Casimir force on a 
static spherical boundary, it is clear that a compensating elastic pressure will 
be built up on the inside. Because of this, most possibilities for measuring 
the surface force that one may conceive of, would not be able to detect the 
electrostrictive contribution at all. A measurement of electrostriction would 
require a measurement of the local pressure in the interior region. 

Consider finally the third term in (1). This is the Abraham term, de- 
tectable under special circumstances at low frequencies |Q, but fluctuating 
out at higher frequencies, especially at optical frequencies. In our case, the 
Abraham term does not contribute. It fluctuates out. 

Summing up so far, we see that it is the Abraham - Minkowski term in 
(1) that gives the most important contribution to the surface force. We will 
be concerned with this term in the following, and will consider both sub- 
classes (1) and (2) listed above. In the next section we highlight the general 
Green-function formalism, and explain the regularization method that we 
use. Section 3 is devoted to a study of the nondispersive force. By using the 
expansion technique mentioned above we derive the expression for F nondls P^ 
to order 1/V 4 , in (27). This force is repulsive. Some earlier results, derived 
by purely analytical means, are corrected. As regards the dispersive force, 
considered in Section 4, we restrict ourselves to giving essentially estimates. 
For high ultraviolet nondimensional cutoff xq the dispersive force is given as 
the sum of expressions (39) and (46). The dispersive force for a compact ball 
is strong, and it is repulsive. It may be as large as about 10 6 — 10 9 times 
the nondispersive force. This fact does not seem to have been universally 
recognized before. 

We put henceforth h — c—1. Electromagnetic Heaviside - Lorentz units 
are employed. 



2 Basic formalism. Regularization 

The general formalism was developed in our earlier paper [1], but for the 
benefit of the reader we will recapitulate some of the essential points. 

The dispersion of the material is accounted for in a simple way, by taking 
the permittivity e(iu) as a function of imaginary frequencies a) to be a step 
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function: e(iu) = n 2 = constant > 1 for uj < ujq, e(iu) = 1 for to > ujq. 
Thus there is one single "absorption" frequency cjq, serving as an ultraviolet 
cutoff. 

When considering the two-point functions, we let the two spatial points 
r and r' be separated in the radial direction. There are two scalar Green 
functions in the problem, Fi(r,r') and Gi(r,r'), defined by 

r,r f <a: F h G t = inkji(nkr < )[hf\nkr > ) - Ap^k^j^nkr^}, (2) 



r, r' > a : Fi, Gi = ik 



3l (k r< ) - B F , G {ka)hf\k r< ) /i, w (At>). (3) 



Here k = \u\, j\ is the spherical Bessel function, and hf^ the spherical Hankel 
function of the first kind. The coefficients A FiG and B FG are related to the 
boundary. 

We shall work with the two-point functions that contain the surface - 
induced contributions only, and thus have to subtract off the contact terms 
corresponding to a uniform medium. It means that the two-point functions 
refer to "disturbed" quantities caused by the boundaries. Specifically, we let 
(i) the contact force on the inner side r = a— be calculated as if the inner 
medium be filling all space; (ii) the contact force on the outer side r = a+ 
be calculated as if the outer medium be filling all space. This means that 
the first term in (2) has to be subtracted off as a contact term. Similarly the 
first term in (3) has to be subtracted off. 

We stress that this method of regularization is not an arbitrary choice. It 
is the only kind of regularization that separates off the volume terms, and it 
is the only method that permits us to write the relation between the surface 
force F and the Casimir energy E as F = —(l/Ana^dE/da. It implies, in 
particular, that both the surface force and the Casimir energy goes to zero 
when the radius goes to infinity, what seems to us a most natural result. Our 
method of regularization is the same as used by Milton jn|, |Tj|, [jl8| ; cf. 



also [14 , |17| . It should be mentioned, however, that this method is different 
from the one advocated recently by Carlson et al. |19[, |F2D[. We shall call the 
modified two-point functions, constructed by the boundary-related parts of 
the Green functions only, the "effective" two-point functions, and write them 
as <> e ff. We introduce the nondimensional cutoff frequency 

x = uj a, (4) 

and perform a complex frequency rotation 

k — > ik = ilu, ka — > iua = ix, (5) 
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whereby the Riccati - Bessel functions si and e\ defined by 

si(x) = J^I u (x), e t (x) = yJ^K v {x) (6) 

correspond to the Wronskian W{si,ei} = — 1. Here v — I + 1/2, l v and K v 
denoting modified Bessel functions. 

The effective electric field products in the limit r' — > r in the interior 
region r < a are 

(E r (r)E r (r')) eff = -— jf — £ -^Z(Z + 1)A G (^ (rur), (7) 

(E ± (r)E ± (r')) e// = -Lj r ^ £ ^-tl{A F (x) S f(^f) - A G (x)[ Sl (nxf-)r}, 

(8) 

with r = r/a, n = n(ix), prime meaning different at ion with respect to the 
whole argument. We need the expressions for the frequency- rotated coeffi- 
cients: 

. , , ei(nx)eii(x) - nei(x)eii(nx) . . 

A f{X) = —, : : ^ z r, (9) 

si[nx)ei![x) — nei[x)si![nx) 

nei(nx)eif(x) - ei(x)eif(nx) 

A g{X) = -, : t\ T\ 7 7 ( 10 ) 

nsi{nx)ei/(x) — ei{x)si/{nx) 

(it may be noted that Ap^ix) = (— Af,g{x)) ■ Similarly in the exterior 
region r > a we have 

1 f x ° dr 2/ -4- 1 
(E r (r)E r (r')) eff = -~ 4 jf - £ -±-Z(Z + l)B G (x)e 2 (xr), (11) 



(£ ± (r)£ ± (r')> e// = -b / ^ E ^{5 F (x)e?(xf) - ^)M*f)] 2 }, 



/ xdx / 

nr 2 a 2 J ^ An 

(12) 

with 

„ , . Si(nx)sif(x) — nsi(x)sif(nx) , . 

B F (x) = ' , 13 

si(nx)ei/{x) — nei{x)si/(nx) 
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B G (x) 



nsi(nx)si/(x) — si(x)si/(nx) 
nsi(nx)eif(x) — ei(x)sif(nx) 



(14) 



The surface force density, as calculated by use of the Maxwell stress tensor, 
is F = F int + F ext , where the interior force is 



int 



- / nxdx 

J{J 1=1 



2l + 1 A M 

4tt ' 



ln( 



sif(nx) 
siinx) 



with 



Af 1 * = [A F (x) + A G {x)\si{nx)sii{nx) 
Similarly the exterior force is 



Fpx. 



ext 



2na 4 



xo 



^21 + 1 



i=i 



47T 



A ext 



ln( 



-eif(x) 

ei(x) 



with 



(15) 



(16) 



(17) 



Af* = [B F (x) + S G (a;)]e,(ar)e,/(z). (18) 

These expressions assume the step-function dispersion relation for n(ix), as 
explained above. The value of the constant n for x < Xq is however at this 
stage arbitrary; (n — 1) is not necessarily small. 

The natural question is now: are the force expressions (15) and (17) finite? 
The answer is no. This can be seen by means of the uniform asymptotic series 
to be employed below, or by direct machine calculation. Thus regularization 
is called for. The regularization procedure is quite straightforward, at least 
for the case of dilute media, the case to which we now turn. 



3 Nondispersive force 

We assume henceforth (n— 1) to be a small quantity, and expand the formal- 
ism, by means of Maple, up to order (n — l) 2 . We make use of the uniform 
asymptotic expansions for the Riccati-Bessel functions, as explained in |l|. 
As mentioned above, we insert these expansions up to order 1/u 4 . 

As we mentioned in Sect.l, the nondispersive case can be treated by 
setting xo = oo. We introduce the symbols 

z = x/v, t(z) = (l + z 2 )- 1/2 , (19) 
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and calculate in the interior region the integrand of Eq.(15), using the Debye 
expansion. Taking the integral from z = to z = oo we obtain 



dznu 3 zA} nt 



ln( 



(n — l)iri> 2 



+ 



bin - lYitv 



2 r 



32 



1 - 



592 1 229 1 



+ 



41 1 
1 ~ 128^2 + 

29504 1 



si'jnx) " 
si(nx) 
711 1 

16384^ 

3771 1 " 



525 iiv 640 v 2 75075 7rz/ 3 81920 z/ 4 



(20) 



This is to be inserted into the expression for the interior nondispersive surface 
force density F™ t ndtsp which, according to (15), can be written as 



-inondisp 
int 



An 2 a 



j oo 



1=1 



dznu 3 zAi nt 



ln( 



sif(nx) 
sAnx) 



(21) 



In the exterior region we calculate similarly 



[n 



2 r 



dzu 3 z Af xt 



41 1 



-eif(x) 
ei(x) 
711 1 1 



ln( 



128 v 2 16384 1/ 4 



3(n - Vf-Kv 
32~ 



2 r 



1 - 



592 1 



33 1 



315 ttu 128 v 2 
29504 1 639 1 



45045 iris 3 16384 z/ 4 



(22) 



and the exterior nondispersive surface force density F ( 



nondisp 
ext 



IS 



rpnondisp 
ext 



4vr 2 



j oo 
"2^4 



1=1 



ei{x) 



(23) 



When constructing the total nondispersive surface force density as F nondls P 

^nor, 
int 



pnondisp _|_ pnondisp ^ wg gee from (20) and (22) that the terms of order (n — 



1) compensate each other. The surface force accordingly starts with terms 
containing (n — l) 2 . 

(The application of the general formalism developed in QXJ] to dilute balls 
contains a calculational error; Eq. (94) should have been multiplied by a factor 
1/2. Therefore the final conclusion of Sect. 4 in |TJ is incorrect. We thank 
K. A. Milton for informing us about this point.) 
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Adding (21) and (23) we now get 



pnandisp = (» ~ g 2 _ J^TI (M) 

64tto 4 ^ v 128 16384 v 2> y ' 

i=i 

Here the divergence mentioned above is made explicit; regularization is needed 
in the first two terms. We shall use the Riemann zeta function regularization. 
This method has in recent years become a standard tool in all areas of quan- 
tum field theory, and is found to be in agreement with what one can derive 
by other means, such as the dimensional regularization, cutoff regularization, 
etc. (cf, for instance, the monograph of Elizalde et al. |2l]]). Zeta function 
regularization has been applied in the context of spherical dielectrics in Jl|] 
and @. An interesting contribution to this technique is further provided by 
the recent paper of Lambiase et al. [P2]| . Here, no formal divergences appear; 
an extra parameter s is introduced whereby all expressions are retained finite. 
At the end, an analytical continuation of s is performed. The final results 
are in agreement with what one finds by a straightforward application of the 
Riemann zeta function technique. 

To effectuate this kind of regularization we need in practice only the 
formula 



X> s = (2- s - 1)CH0, (25) 



1=0 



according to which 



oo 1 oo 

5> 2 = 4 £-° = -l- (26) 



4 

i=i i=i 



The third sum in (24) is finite and is known exactly: YlHi v ~ 2 = ( 7I " 2 /2) — 4. 
From (24) we then get, when writing the force as F nond%sp \i/ u i to emphasize 
that it is based upon a uniform asymptotic expansion to order 1/z/ 4 : 

pnondisp\ ^ _ ( n ~ jj r ^29 927 

11/1/4 ~ 4vra 4 [ 65536 524288 J ' 1 ' 

It is of interest to compare this expression with the expression F nondlsp \i/ u 2 
following from a uniform asymptotic expansion to order 1/z/ 2 only: 

F nondispl = -J- — . (28) 

17 4?ra 4 2048 V ; 

The reason for this comparison is the following: the combination of terms 
occurring in the surface force density (24) involves even powers of 1/z/ only. 
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The structure of terms is such that the 0(l/v 4 ) force is somewhat larger than 
the correct force, whereas the 0(1/ v 2 ) force is somewhat smaller. From (27) 
and (28) we are therefore able to give quite accurate bounds for the "exact" 
nondispersive surface force density F nondls P ; 

pnondisp^ j 2 <i pnondisp ^ pnondisp^ ^ (29) 

(The term "exact" here refers to the uniform asymptotic expansion only; 
we are of course concerned with a dilute- medium theory of order (n — l) 2 .) 
Numerically, 

(n - I) 2 (n - I) 2 

' -0.0161 < F nondisp < ^ — -^0.0194. (30) 



\;;<r iira 4 

The nondispersive part of the surface force is thus repulsive. The quantum 
mechanical result is in this sense in agreement with the prediction of the 
Abraham-Minkowski term in (1). 

We note that some care has to be taken to include a sufficient number of 
terms of the uniform asymptotic expansion when constructing the expression 
for F nondisp . To illustrate this, let us see what becomes the result from 
including only the first term in (24) (this corresponds to including only the 
zeroth order of the uniform asymptotic expansion). Making use of the first 
regularizing equation in (26) we obtain 



(V? _ i V i 

= --rJ-si- (31) 



Since 1/64 = 0.0156 we see that the magnitude of (31) is not very different 
from the second order approximation in (30), but the sign is wrong. We 
conclude that it is necessary to go to at least the second order in the uniform 
asymptotic expansion to get the sign right. The expression (31) is actually 
in precise agreement with the nondispersive force as calculated by Milton 
||15|| , [|l~6|l , and Milton and Ng [Q. The reason why these authors got a 
discrepancy with (30) is simply that they included an insufficient number 
of terms . In ]17|] and the following accurate value is given, with the 
assistance of numerical methods, 

(n - I) 2 

pnondisp = \ ' .01907 (32) 

The fourth order result given in (27) is about 1.7 per cent high. 
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4 Dispersive force 



This case is more difficult to handle than the nondispersive case, and we will 
essentially be able to give only order-of-magnitude estimates. We restrict 
henceforth the uniform asymptotic expansion to the second order. Inclusion 
of higher order terms would here be pointless. The ball is assumed to be 
dilute, as before. 

We now have to keep the nondimensional ultraviolet cutoff Xq = ujqO, as a 
finite quantity in thr formalism. One may ask: what are the physically rea- 
sonable magnitudes of x$ ? We note that xq is here an "external" parameter 
that has to be inserted into the formalism by hand. Its value is given by the 
cutoff frequency ujq , which is determined by the molecular structure of the 
medium, and by the radius a. We will take uq to lie in the ultraviolet re- 
gion. It seems reasonable, by comparison with the Lorentz dispersion model 
(cf. [fffil ) for instance), to put u = 3 x 10 16 s _1 . Thus in dimensional units, 
assuming a small ball with radius a = lfim, we get xq = lo^ajc = 100. Balls 
with somewhat larger radii can also be actual. Accordingly, 

100 <x < 1000 

seems to be a reasonable range for xq. 

Another question is: what is the upper physical limit I — Iq to be taken 
in the summations? This point becomes accentuated by the circumstance 
that some of the sums are formally divergent. Obviously, this behaviour is a 
spurious effect. Our dispersion relation, as explained in the beginning of Sect. 
2, implies that photons having frequencies uo higher than uq do not "see" the 
medium at all. If a photon of limiting frequency u just touches the surface of 
the sphere, its angular momentum is equal to u a , i.e., equal to x . We thus 
expect, on physical grounds, that l is of the same order as x . This kind 
of argument has repeatedly been used in Casimir calculations ||, 0, [23 . 



It means that one confines oneself to giving order-of-magnitude estimates 
rather than exact numbers. The physical force, of course, is a definite number 
in a given case, but to evaluate it one needs detailed information about the 
dispersive structure of the medium. That lies outside the scope of the present 
paper. 

Adding (15) and (17), making use of the second order approximation, we 
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get 



(n-l) 2 ~ - r*o/* 



47r 2 a 



4 



1=1 J ° 



i r x o/ u i 
+ / zV2 W( 3 - 42 ^ + lfe4 + 2z& ) dz + ( 33 ) 

V 

Here the first and the second integral stem respectively from the zeroth, and 
the second, order terms in the uniform asymptotic expansion. We write F 
as a sum of two terms, one zeroth order term A and one second order term 
B, and consider B first. This term is more easy to handle than the term A, 
since the sums over I in B are convergent and can be calculated by means of 
the Euler-Maclaurin formula 101. It is convenient to write B in the form 



(n - l) 2 1 ^ r xo 2 3v 9 - 42xV + 18x 4 z/ 5 - 2x 6 z/ 3 
~ 4vr 2 a 4 4 J % (x 2 + v 2 f : 

and to use the Euler-Maclaurin formula in the following version: 

L 

dlf(l) 



£/(o = £/co + / 
i=i i=i Jl+i 



+±[/(oo) + /(L + 1)] + l[/'(oo) - f (L + 1)] + ... (35) 

Here L is an auxiliary integer, the use of which may improve the accuracy of 
the series. Usually it suffices to let the magnitude of L be moderate (L < 10; 
cf. [0] for instance, where we made an analogous calculation). We write the 
expression (34) as B = Yli=i /(0> ano - obtain for the first term in (35) to a 
good accuracy 



L 

£ 

i=i 



f(l) = ^^J j J™ z 2 t 12 (z)(3 - A2z 2 + lSz 4 - 2z 6 )dz 



in ~ 1)2 " h -L. (36) 



4vra 4 2048 



Here we could set the upper limit of the integral equal to infinity, since 
the main contribution comes from z close to the lower limit. The second 
term in (35) yields after some algebra, when introducing the abbreviation 
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w = L + 3/2, 

dlf(l) 



L+l 



n — l) 2 x 1 w s x 



4vr 2 a 4 L 8 64 (w 2 + x 2 ) A 
\ 1 w 6 x 105 w 4 x 



128 (u> 2 + :r 2 ) 3 512 (u> 2 + X 2 ) 2 
325 w 2 x 315 rr 125 



H u> arctan wn\ 

■ w 1024 J 

wir), x — > oo. (37) 



1024 w 2 + a;g ' 1024" ~ w 1024 
(ra - l) 2 x 65 



47r 2 a 4 v 8 2048 

The last term between square parentheses in (37) was evaluated putting 
x = oo. Similarly we put x = oo when evaluating the last nonvanishing 
terms in (35): 

+ + (38 ) 

Altogether, for large cutoffs, we get the L-independent result 

„ (n — l) 2 , Xn 65 . . 

B = w<-f + aoS**' (39) 

We extract the nondispersive part of this expression: 

,. . x (n — l) 2 65 , in . 

B [nondispersive part) = — — j — 2048 ' 

This is compared with the second order contribution to the nondispersive 
force obtained previously in Eq.(24): 

F non d is P(2nd Qrder) = (^Zi0! ( __^) y u = (n-lf 65 

v ; 64vra 4 v 128 4vra 4 2048 v ; 

i=i 

The expressions (40) and (41) are seen to be equal. This brings us to the 
following important conclusion: If we perform a dispersive calculation of the 
surface force, to second order, we find a term that is independent of the 
cutoff. This term is precisely equal to the second order contribution within 
a calculation that is nondispersive from the outset. In the second case zeta- 
function regularization is being employed, while in the first case it is not. 
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Consider next the zeroth order term A in (33). Performing the integration 
over z we get 

n-1) 2 ^. „ r° /u 



On — i Y 2 r Xa/u 



1=1 

(n-1) 2 ^ 1 z/ 5 x 1 v 3 x v 2 x 

= - — tt-t- > — ~ f — « ^777 H o o H arctan — . 42 

47r 2 a 4 ^ L 4(z/ 2 + x 2 ) 2 8z/ 2 + a;g 8 v J V ; 

We will only be able to give an order-of-magnitude analysis of the physically 
meaningful part of this divergent expression. Note first, however, the follow- 
ing point: take the nondispersive limit, Xq = oo. Then only the last term in 
(42) survives and we get, again employing zeta function regularization, 

(n — 1) 2 1 ^ (n — l) 2 1 
A(nondispersive part) = —— —— v 2 = — — -. (43) 

i=i 

This is precisely the same as the expression (31), and serves as a corollary of 
our dispersive calculation. 

Now return to the dispersive case, and estimate the contribution to the 
sum in (42) from low and moderate values of /, such that v < x . When 
i/<xo, we see that the two first terms in (42) are negligible in comparison 
to the third term. Moreover, when v = x the two first terms compensate 
each other. Consequently we may write 

a\ v < xo ~ Yl j arctan -■ ( 44 ) 

u<xo 

Here the variation in the factor arctan(xo/^) is moderate; it decreases from 
7r/2 for small v , to 7r/4 when xq = v . For estimate purposes we may replace 
arctan(xo/^) by a constant C, where C lies in the interval 7r/4 < C < n/2. 
The sum in (44) then becomes simple to evaluate, since 

*°-V2 1111 

" 2 = o x o + 2 X ° + 6 X ° ~ 4' (45) 

i=i 



and so we get 



{n-lfC A , 1 o 1 1 . 



It is here natural to argue such as we did above, in connection with (40) 
and (41). We may thus associate the x - independent part of (46) with the 
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nondispersive surface force, calculated according to the zeroth order uniform 
asymptotic expansion. Actually, if we equate the x - independent part 
(arising from the last term in (46)) to the expression (31), we obtain C = 
7r/2. This is precisely the value of C corresponding to the small values 
of v , according to our discussion above. We may thus conclude that the 
nondispersive part of the force, at least as far as the leading term is concerned, 
is associated with the low angular momenta I . 

When x is included, (46) predicts under usual physical circumstances 
a strong, repulsive force. For instance, if we take C = n/2 and Xq = 100, 
we find that (46) is about 10 6 as large as the nondispersive second order 
(or fourth order) force given in (30). This large difference can be expected 
to be quite important, in practice. An eclatant example to think about is 



the sonoluminescence effect [[14]]. We will briefly return to this point in the 
Conclusion section. 

The result (46) was derived under the assumption of low and moderate 
nondimensional frequencies, v < xq. Is there an appreciable contribution to 
A from higher frequencies also? We think that this is not so, the reason being, 
as explained above, that photons with u> > u do not "see" the medium. It 
is physically speaking most safe to truncate the I summation at v ~ x 
and adopt (46) as a reasonable order-of-magnitude estimate of the leading 
dispersive term. The total surface force density F is thus estimated adding 
the expressions (46) and (39). 

Finally we mention the following alternative method of handling the di- 
vergent sum A, which may appear mathematically convenient. The method 
consists in adding and subtracting the weakly divergent sum representing A 
in the frequency region v > xq. Starting from (42) we see that this amounts 
to writing 

(n-1) 2 ^ 1 v 5 x 1 z/ 3 x v 2 x x\ 

A = 7T-r- > — tt, ^ttt H n o H arctan + 

4vr 2 a 4 j^ 1 A(v 2 + xl) 2 8 v 2 + x 2 Q 8 v 3z/ J 

(n-lf xlyl (47) 

where we have chosen to let the sums over I in the extra terms run from 
1 to infinity. The advantage of writing A in this way is that the first of 
the two sums in (47) is finite, and that we have explicit control over the 
divergence in the second sum. Let the first sum be denoted by A (finite). 
Making a numerical calculation of this expression at the lower and upper 
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limits xo = 100 and xo 

A{ finite) : 



1000, we get approximately 
(n- I) 2 



An 2 a 4 



-1.7 x 10 6 , 
-2.4 x 10 9 , 



x 

X 



100 
1000 



(4£ 



A graphical representation of A(finite) versus Xq shows that the variation 
is roughly linear in the actual frequency region. Exploiting this fact, and 
truncating the second sum in (47) at I = xq, we arrive at the following 
estimate for A: 



A 



{n - If 
An 2 a 4 



x 



3 x o 



^£ I -(266.5 

l=i 



100 



264.8U0 6 



(49) 



valid for 100 < xq < 1000. It is seen, in fact, that the structure of (49) 
is not very different from the structure of our previous estimate (46). The 
advantage of (49) in comparison with (46), of course, is that there appears 
no extra constant C to be estimated separately. 



5 Conclusion and final remarks 

The physical model considered in this paper is that of a compact dilute spher- 
ical ball, surrounded by a vacuum. The dispersion relation of the material 
is such that n(iu) = n = const when the imaginary frequency Cj is lower 
than the ultraviolet cutoff frequency Uq whereas n(iCj) = 1 for higher Cj . 
This is a very simple dispersion relation: we know on general grounds that 
n(iCj), being closely related to a so-called generalized susceptibility, is real on 
the positive imaginary frequency axis and decreases monotonically towards 
1 when Cj — > oo (cf., for instance, f2"5|). However, this is not a serious restric- 
tion here as we may well imagine that n{iuj) instead is endowed with a small 
negative slope for Cj < u . 

The nondimensional cutoff is Xq = ojQa, a being the radius of the ball. 
Under typical physical conditions, if the radius lies in the interval lum < 
a < lOfim, we expect that 100 < xq < 1000. 

Our basic Green function formalism (Sect. 2) is Taylor expanded in 
(n — 1) up to the second order, and is expanded in 1/u up to the fourth 
order [v = I + 1/2). The Maple V program here proves to be an effective 
tool. 

The surface force density F is constructed by taking the difference be- 
tween the radial Maxwell stress tensor components on the two sides of the 
boundary r = a . Alternatively, F can be found as the integral of the radial 
component of the volume force density f across the boundary. The general 
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expression for f is given in (1), and the meaning of the various terms is dis- 
cussed. Only the first, and under usual conditions most important, term (the 
Abraham - Minkowski term) is further considered in this paper 

The nondispersive force is calculated in Sect. 3. Formally, this case is 
treated by setting x = oo, and regularizing the divergent sums using the 
Riemann zeta function. Working to the order 1/z/ 4 in the uniform asymptotic 
expansion, the nondispersive force is given by (27) and (30). A slightly more 
accurate value is given by (32). This force is repulsive. 

We emphasize that F nondls P is a physical force; the expression is not merely 
a formal outcome of a more or less arbitrary regularization procedure. It is 
here of interest to note that, in addition to the Green function approach 
above, there has recently been developed also two other alternative proce- 
dures which lead to results supporting our statement: 

First, one may consider the problem from a microscopic point of view and 
calculate the mutual van der Waals energy for the molecules in the compact 



sphere. This has been done by Milton and Ng pq| . The van der Waals energy 
for two molecules is V = —B/r 7 , with B a positive constant. Integrating this 
over the sphere one finds a divergent expression which has to be regularized, 
conveniently by means of dimensional regularization. The finite remaining 
part of the energy becomes 

E^-^l^ (50) 
vra 1536 v ; 

Associating the van der Waals energy E with the surface force F via the 
relation F = —(l/Aira 2 )dE/da, we find the same numerical result for F = 

Second, one can calculate the mutual energy between the molecules in 
the sphere by employing quantum mechanical perturbation theory. This has 
been done by Barton |27], for a dilute sphere. He used an exponential cutoff 



for the wave numbers. Again, the formal expression for E was found to 
diverge, but after omission of the cutoff dependent terms, the result was 
found to be precisely in agreement with the expression (50). 

These two recent alternative developments are welcome, since they signify 
that a consensus between various approaches to the difficult Casimir surface 
force problem is finally in sight. 

We emphasize that the repulsive force F nondls P^ although being a physical 
force and thus measurable in principle, is not the total surface force, not 
even for a nondispersive medium. The total force is necessarily attractive, 
since it results from the attractive van der Waals force between molecules. 
(This point was apparently regnogized very early, by Davies in his classic 
paper [EH].) At first sight it may perhaps be surprising that we operate with 
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a measurable force component which is only a part of the total force. This 
kind of situation is however known also from other areas in electrodynamics, 
and it by be instructive in the present context to discuss briefly the following 
two examples: 

1) Consider the situation where a high intensity laser beam falls vertically 
from above on a liquid (water) surface. The surface bulges out, corresponding 
to an outward-directed force. The elevation is slight, of the order of 1 \im, 
but is clearly detectable. This is the famous Ashkin- Dziedzic experiment 



29| . The elevation can be described completely, using only the Abraham- 
Minkowski force term, i.e., the first term in (1). 

The Abraham-Minkowski force is however not the only force acting on the 
water surface. There is an additional force, arising from electrostriction, and 
this force is stronger than the first-mentioned one such that the total force on 
the surface is actually compressive. Electrostriction does not contribute to 
the elevating force, as mentioned earlier, since electrostriction merely com- 
presses the column of liquid. The theory of the Ashkin-Dziedcic experiment 
has been given in [30| and [ffl. 

2) A second, related, situation is the one discussed in many textbooks ( 
cf. ||, for instance), namely two parallel condenser plates partly immersed in 
a dielectric liquid. When a strong static electric field is applied between the 
plates, the liquid is observed to rise slightly in the interior region. The situa- 
tion is more involved than one might be inclined to think. Again, the observ- 
able elevation can be descibed entirely with use of the Abraham-Minkowski 
force only. There is in addition the compressive electrostriction force acting, 
so that the total surface force density is compressive. The electrostriction 
force has to be the stronger member of the two forces, since otherwise the 
column of liquid would not be kept together as a whole. The mathematical 
reason why the electrostriction force is the stronger one, is that the relative 
permittivity is larger than 1. This is a thermodynamical result. Seen from 
this perspective, the deepest reason why the column of liquid is kept together, 
is neither hydrodynamics nor electrodynamics, but rather thermodynamics. 
The theory of this situation is discussed, among other places, in M and |JT] . 

The lesson to be learned from the above two examples is fairly obvious: 
it is physically quite legitimate to consider a part of the complete surface 
force, and discuss its experimental consequences, without drawing the com- 
plete force into consideration at all. This is precisely what we do in principle 
when we discuss the Casimir force F™"^. We are not aware of any con- 
crete operational argument so far to really measure F nond%s v experimentally, 
however. 

Consider next the dispersive force. It was calculated in Sect. 4, to second 
order in the asymptotic expansion. The force is written as a sum of two 
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terms, F = A + B, where the second order term B is given in (39). This 
term, although strictly speaking derived for high cutoffs only, xq — > oo, is 
expected to be fairly accurate for all values of x that are actual in practice. 
The term B has the same basic structure as the surface force calculated 
earlier on a spherical shell 0: there is an attractive part linear in Xq, which 
is due to dispersion, and in addition a repulsive part which is dispersion 
independent. Under usual physical conditions the first of these two parts is 
the stronger one. 

As for the leading term A, of zeroth order, we have only been able to 
give an order-of-magnitude analysis, based upon physical arguments allowing 
us to truncate the divergent sums over I at an upper limit approximately 
equal to xq. We gave two estimates for A, in equations (46) and (49). The 
structures of these two expressions are seen to be roughly similar. Under 
usual circumstances also the force A is a repulsive one. 

To get an idea about the magnitudes involved, let us limit ourselves to 
the dominant first term in (46), putting for definiteness C = n/2. Then 

pdisp „ (n-lf xl = {n-lf ul 

4vra 4 48 4vra 48' K ' 

This force may easily be about 10 6 as strong as the nondispersive force (27). 
The force (50) corresponds to the Casimir energy (integrating while keeping 
ujq constant) 

E dis P „ _ {n-l) 2 u>l 2 

96 v ; 

This energy is always negative, and has its maximum when the radius is at 
minimum. Moreover it is independent of the substitution n — > 1/n, when n 
is close to unity. One would therefore think that the formula is equally well 
applicable to the case when there is a vacuum bubble immersed in a dilute 
dielectric medium. 

One may wonder whether it is possible to apply the present kind of the- 
ory to the phenomenon of sonoluminescence (an instructive review of this 



phenomenon has recently been given by Cheeke [32 1. It is well known that 
the energies involved in the nondispersive Casimir effect are far too small to 
account for the sonoluminescence effect. Can the very large magnitude of 
(52) "save" the situation? We see that by taking (n — l) 2 ~ 0.1, x ~ 1000 
(corresponding to a ~ 10/im), we obtain from (52) E d%sp ~ —2 x 10 4 eV. The 
magnitude of this energy is probably acceptable, but the sign is apparently 
wrong fl~l, |4f. At least the sign is in conflict with energy emission during 



the collapse. Consider, however, the following argument: an experimental 
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fact is that the photon flash does not occur during the whole collapse period 
but instead precisely at the end of that period, when a = a m i n ~ 1/j/m. Imag- 
ine now that the energy of the burst of photons were really taken from the 
Casimir energy at this particular instant. Consequently there would occur 
an energy imbalance in the system. The balance of Casimir energy would be 
restored if the bubble were increasing its radius out to a ~ lOfim. 

Such a movement actually happens. Evidently, we do not "explain" the 
sonoluminescence effect by this kind of argument. Our aim is more modest, 
viz. to suggest that sonoluminescence is perhaps not totally unrelated to the 
Casimir effect after all. 

The last point that we shall dwell on in this paper, is the influence from 
absorption in the medium. This is an important point since it shows the 
physical limitations of our adopted material model, and the directions for 
further research. At first sight, it would seem as if absorption does not cre- 
ate any problems at all in the theory: we are after all adopting a model in 
which the quantity n(iuj) is real and shows a simple though physically ac- 
ceptable behaviour along the positive imaginary frequency axis. Absorption 
should accordingly be taken into account implicitly in the theory, through 
the Kramers-Kronig relations, which lead to a complex behaviour of n{uo) 
along the real frequency axis. The problem waiting for a complete solution 
lies at another place, namely in our adoption of formulas (7) and (8) for 
the effective products. These are the same formulas as for a nondispersive 
medium. A more complete theory would have to take into account the cou- 
pling of the molecules to a heat reservoir, with a corresponding quantization 
of the source currents. This would in turn imply the presence of Langevin 
terms in the effective products. Interesting developments of the theory of 
absorptive media, implying a quantization of the complete system (medium 
plus reservoir) have been given in P5[ , |34|| , and most recently, in |35| and 



|36| . We only point to this problem here, without going into any detail. 
For the convenience of the reader, we summarize in Appendix B the basic 
assumptions leading to our formulas for the effective products. 
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A Nondispersive force expanded to the fifth 
order in (n — 1) 

Our calculation of the nondispersive force in Sect. 3 was based on a dilute-ball 
expansion up to order (n — l) 2 , and on an asymptotic expansion up to order 
1/u 4 . A characteristic property was found to be that terms of order (n — 1) 
in F™ i ndlsp and F™°/ ldlsp compensated each other, so that the leading term in 
pnondisp = F ^ndis P + pnondisp wouU gQ as ^ _ rf ^0) and (22). One 

may ask: what happens if one makes a higher order expansion in (n — 1), will 
then terms containing odd powers in (n — 1) continue to compensate each 
other? We shall investigate this point, giving the expansions for F^ dlsp 
and F™™ dlsp up to order (n — l) 5 . However, we will now have to truncate 
the uniform asymptotic expansion at order 1/u 2 . (A combined set of higher 
order expansions in (n — 1), as well as 1/u , turned out to be beyond the 
capacity of our Maple program.) 

We abstain here from giving the explicit expressions for F™ t ndlsp and 
Fext d%SP : but write down their sum: 

pnondisp = ( n ~ 1 ) 2 r 33 _ / _ jW 99 , 11 ^1 

4vra 4 L 2084 V A 4096 210 tt 

1N2 . 1213 11 1, . ^3. 785 29 l Wco , 

+(n - 1) ( 1 ) - (n - 1Y( 1 ) (53) 

V ; V 65536 105 tt ; v ; v 131072 216 n ny ' 

The first term in this expression agrees with our previous second order result 

(28). We see that odd powers in (n — 1) generally do not compensate each 

other in the surface force. All powers in [n— 1) are present, with the exception 

of the first, linear, term. The expansion given above may be useful also in 

cases when (n — 1) is not so small, perhaps up to (n — 1) ~ 0.3 — 0.4. 

B The effective products derived from statis- 
tical mechanics 

Assuming for definiteness Schwinger's source theory, we relate the electric 
field components Ei(x) to the polarization components Pk{x') via a tensor 

r«fc ) . 

E t (x) = J d 4 x> T tk (x, x') P k (x') (54) 
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Stationarity of the system means that depends on time only through the 
difference r = t — t' '. Causality implies that the integration over if is limited 
to t' < t. From a statistical mechanical viewpoint, r^(x, x') is a generalized 
susceptibility p5 |. 

Introduce the Fourier transform r^r, r', a;) via 

T tk (x, x') = f°° P e-^ T tk (r, r', u) (55) 

The basic equation (54), in combination with the causality principle, are thus 
the necessary conditions for causing the function T ik (r,r' ,u) to satisfy the 
following properties: 

(i) It is a regular, one- valued function in the upper half frequency plane. 
At infinity (\u\ — > oo), it goes to zero. 

(ii) There is no singularity on the real axis (except at the origin, in the 
case of a metal). 

{Hi) The function does not take real values at any finite point in the upper 
half plane except on the imaginary axis. 

Now invoke Kubo 's formula from statistical mechanics: it states that the 
spectral generalized susceptibility is related to the commutator between field 
components as 

poo 

r, fc (r,r» = i / dre^ T ([E^x), E k {x')}) (56) 

(cf. for instance, |p5]). This means that the generalized susceptibility is 
identifiable with the retarded Green function: T ik (x, x') = Gf k (x,x'). For 
t < t' both T ik (x, x') and G^ k (x,x') vanish, the first because of causality, the 
second because of the definition of the retarded Green function. 

Consider now the two-point function (Ei(x)E k (x')) . Its Fourier transform 
{Ei(r, u)E k (r', uj')) can be expressed in terms of the quantity (i?j(r)_E , fe (r')) w 
(the spectral correlation tensor) as 

(^(r^^rV)) = 27f(E i (r)E k (r')) UJ 5(uj + uj') (57) 



The crucial step now is to apply the fluctuation- dissipation theorem |25| . 
according to which 

(^(r)E fc (r')). = ImG^(ry,u) coth(^), (58) 

with P = \/{k B T). At zero temperature, coth(|/5u;) — > sgn{uj). Choosing 
the causal integration contour in the complex frequency plane we obtain, 
after a complex frequency rotation, the effective products as given in Eqs.(7) 
and (8). 
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